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ABSTRACT 



Two-dimensional turbulent wakes in constant adverse and 
favorable pressure gradients were studied experimentally. 

The static and dynamic pressures through the wake were ob- 
tained with calibrated Prandtl-type pitot-static probes, 
and a computer program was written to extract from these 
data the velocity profiles, displacement thickness, momentum 
thickness, and the shape factor. This program further pro- 
vides the turbulent shear coefficient, eddy viscosity, and 
Prandtl mixing length within the wake. 

A wind tunnel was designed specifically for this experi- 
ment and its construction details are included. The tunnel 
test section is unique in that it can be adjusted to generate 
any specified adverse or favorable pressure gradient. 

The results of the experiment indicated that the eddy 
viscosity was a function of the history of the flow rather 
than of the local mean flow field, as is classically assumed. 
The reaction of the turbulent wake to a favorable pressure 
gradient after encountering an adverse pressure gradient was 
significantly different from the reaction of the wake to 
only a favorable pressure gradient. 
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NOMENCLATURE 



A = cross-sectional area 
| | = indicates absolute value 

A^ - turbulent exchange coefficient (Eq. 3, Introduction) 

a = constant in Prandtl’s eddy viscosity (Eq. 5, Intro- 
duction) 

b = width of mixing region (Eq. 5, Introduction) 

3 = constant in test section area distribution 

(Eq. 7, Appendix A) 

= volumetric flow rate through the test section 

Cp = negative of the test section pressure gradient 
(Appendix A) 

1 2 

= turbulent shear coefficient = x/p-pu e 
d - wake generator diameter 

= indicates differentiation with respect to x 

g 

= indicates partial differentiation with respect to x 

d X 

6* = displacement thickness 

, . . . .3u 

e = eddy viscosity = x/^— 

dy 

exp = indicates exponentiation 

F = arbitrary example function 

H = shape factor = 6*/0 

i = indicates i direction (i = 1,2,3) 

00 = large distance in the cross-streamwise direction 

j = indicates j direction (j = 1,2,3) 

L = length of pressure gradient 

Z - Prandtl’s mixing length (Eq. 4, Introduction) 
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logarithm to base e 
laminar coefficient of viscosity 
wake centerline conditions 
static pressure 

1 2 

dynamic pressure = ^pV 
fluid density 

stress tensor in Navier-Stokes equations 

mean flow stress tensor in the time averaged 
Navier-Stokes equations 

indicates time 

indicates the characteristic period used for 
time averaging 

Reynolds shear stress tensor 

momentum thickness 

mean velocity in the x-direction 

fluctuating velocity in the x-direction 

mean velocity in i,j directions (i,j = 1,2,3) 

fluctuating velocity in i,j directions (i,j = 1,2 3) 

mean velocity in the y-direction 

fluctuating velocity in the y-direction 

one dimensional mean flow velocity 

streamwise coordinate 

j coordinate 

cross-streamwise coordinate 

lateral position at which C is being calculated 
(Eq. 5, Appendix B) T 



Subscripts 



e = local external flow conditions 

i,j = directionality parameter (i,j = 1,2,3) 

max = maximum value within the wake 

min = minimum value within the wake 

t = shear related quantity 

x,y = indicates partial differentiation with respect to 
x and y 

1 = test section inlet conditions 

2 = test section exit conditions 

Superscripts 

= indicates a time averaged quantity 
* = designates a fluctuating quantity 
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I . INTRODUCTION 



When tangential stresses are applied to a fluid having 
internal friction, shearing motions are set up in line with 
the stresses and in conformity to the shape of the bound- 
aries. Within this flow field, various kinds of secondary 
motions become possible. Regular ones, when they occur, 
are readily accounted for in terms of pressure gradients 
resulting from the curvature of the main flow, usually 
imposed by the shape of the boundaries. Irregular ones, 
called turbulence, are by far the more common, and their 
direct cause is less obvious. Their occurrence does not 
depend on the shpae of the boundaries, but like all secon- 
dary flows they must depend on a generating mechanism within 
the flow itself. According to present evidence the initial 
onset of turbulence occurs suddenly by a breakdown of the 
laminar flow in localized regions. Turbulence is convected 
downstream in the manner of any other fluid property and 
washed away from its point of origin. Repeated breakdowns 
are therefore generally required to maintain a continuous 
supply of turbulence, and instability of laminar flow is 
an essential part of this process. 

In treating turbulence mathematically it is generally 
assumed that the motion can be separated into a mean flow, 
whose components are u^ (Cartesian tensor notation), and a 
superimposed fluctuating motion whose components are u^, the 
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mean values of which are zero. For the flow of an incom- 
pressible Newtonian fluid, this decomposition can be applied 
to the dependent variables of the Navier-Stokes equations 
and time-averaged, the latter defined as 

t+T 

F(x,y) = is / F(x,y,t)dt (1) 

t 



where T is a period of time large compared with the char- 
acteristic time of the turbulent oscillations. The resulting 
time-smoothed Navier-Stokes equations are 



9u^ 

^atT 



a 

ax . 



(a 



ij 



p u i u j 



- pu i u j } 



( 2 ) 



These equations have the same form as the original 
Navier-Stokes equations if o^. is replaced by - pu^u j . 

Thus the equations for the mean flow are the same as the 
original equations of motion, except that there are addi- 
tional virtual stresses, = -pu'^uj , which represent the 

mean rate of transfer of momentum across a surface due to 
the velocity fluctuations. These virtual stresses were 
first introduced by Osborne Reynolds and are known by his 
name [Ref. 4], 

The quantities u| and uj are the velocity fluctuations 

in the i- and j- directions at one point in space and u’u! 

i J 

is a measure of the correlation between the fluctuations. 
Correlations between fluctuations at two points in space 
are also measurable and give the size and orientation of 
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the turbulent eddies. In trying to solve analytically for 
the second-order correlation tensor, the solution is given 
in terms of third-order correlations. Similarly, the 
solution for the third-order correlation is given in terms 
of the fourth-order correlation, and so on. Hence the 
system cannot be closed and it is impossible to solve 
exactly for the second-order correlations [Ref. 1]. 

The various methods for the calculation of the turbulent 
flow field developed to date are based on empirical hypoth- 
eses. These approaches endeavor to establish a relationship 
between the Reynolds stresses produced by the fluctuating 
velocities and the mean values of the velocity components 
and their spatial derivatives. By this means, the virtual 
stresses are given an analytical form which closes the 
system of equations and constitutes the starting point for 
the calculation of the mean flow. 

The first hypothesis relating the Reynolds stress to 
mean flow properties was advanced by Boussinesq and intro- 
duced a turbulent mixing coefficient, , such that the 
turbulent shear stress could be written 



T = 




(3) 



This coefficient corresponds to the viscosity, p, in laminar 
flow and is therefore often called the "virtual” viscosity. 
Equation (3) has the great disadvantage that the virtual 
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viscosity , A , is not a property of the fluid, as y is, but 
depends on the structure of the flow field. 

In 1925 Prandtl proposed his mixing length hypothesis 
in which the turbulent stresses are modelled as follows: 



The mixing length, £, is that distance in the transverse 
direction which must be covered by an aggregation of fluid 
particles, travelling with its original mean velocity, in 
order to make the difference between its velocity and the 
velocity in the new lamina equal to the mean transverse 
fluctuations in the turbulent flow. The mixing length is 
still not a property of the fluid, but it is, in theory, 
a purely local function and it is far simpler to make 
plausible assumptions concerning the mixing length, &, than 
for the virtual viscosity, A . This constitutes the essen- 
tial superiority of Prandtl* s equation over that due to 
Boussinesq. Even so, Equation (*0 is unsatisfactory in that 
the apparent kinematic viscosity, e = -| vanishes at 

points where is equal to zero, that is at points of 

maximum and minimum velocity, which is contradictory to 
experimental evidence. A simpler expression for e, also 
proposed by Prandtl which overcomes this problem is 







£ 



= ab ( u 



max 




( 5 ) 
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and is valid only for free turbulence, where a is an experi- 
mentally determined, dimensionless constant and b is the width 
of the mixing zone. Thus e, for all the classical models of 
Reynolds stresses, remains constant over the whole width of 
the wake and is a function of the local mean flow. 

Several other similar hypotheses have been proposed as 
analytical expressions for the Reynolds stresses which will 
only be mentioned here. They include: Taylor’s vorticity 
transport theory. Von Karman ’ s similarity hypothesis, and 
the universal velocity distribution law, each of which is 
linked to the mixing length concept due to Prandtl. Each 
of the aforementioned theories is valid only for free turbu- 
lent flow. While none of the theories adequately models the 
turbulent process in general, they are useful tools in the 
integration of the equations of motion in special cases. 

In this investigation it was desired to ascertain, at 
least qualitatively, whether the virtual viscosity could bf 
considered to be a function of the local mean flow propert es 
in the case of two-dimensional turbulent wakes in constant 
favorable and adverse pressure gradients. A wind tunnel 
was designed and fabricated for this study which incorporated 
a variable-geometry test section that could be set to provide 
any specified pressure gradient. While this study was carried 
out using pitot-static probes as data sensors, it is antici- 
pated that the flexibility and utility of the wind tunnel 
test section will allow a very detailed analysis of the 
turbulent wake using a hot-wire anemometer. This will also 



13 



permit the assessment of more modern eddy viscosity models 
which involve hypotheses with regard to the dependence of 
the shear on the turbulent kinetic energy. 
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I I . THE WIND TUNNEL 



A. INLET 

The inlet bellmouth shown in Fig. 1 was constructed 
from eight-inch diameter aluminum pipe halves formed into a 
two-feet square inlet section. The corner seams, which were 
exposed to the inlet air flow, were filled with epoxy putty, 
hand formed, and sanded into a smooth fillet. This was done 
to smoothly transition the flow around these corners into 
parallel flow. 

Immediately downstream of the inlet bellmouth was the 
plenum section of the tunnel. It consisted of four two-feet 
square by one foot long steel sections bolted together. 

Cork gaskets were inserted between the sections to prevent 
the influx of air at the joints and allow for smooth air 
flow over these junctions. Finally, two copper screens of 
#16 mesh between the first two and last two sections reduc d 
the ambient turbulence level in the flow prior to its 
entering the nozzle. 

Initial operation of the tunnel indicated the presence 
of swirling in the air flow through the test section. This 
swirling was eliminated by the insertion of a four-inch 
thick, two-feet square section of quarter-inch honeycomb 
just prior to the first copper screen. 
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B. NOZZLES 



Two different nozzles were constructed and used during 
the course of this investigation. One nozzle, used with the 
favorable pressure gradient test section, was constructed 
from half-inch Plexiglas. The nozzle, shown in Fig. 2, is 
two feet long, has an exit area of 168 square inches, (24 
inches high by 7 inches wide), and a contraction ratio of 
3.42. The top and bottom sides of the nozzle are parallel 
while the vertical sides contract to form the nozzle. The 
wake generator was attached to the nozzle in its exit plane 
on the axial centerline of the vertical walls. Static 
pressure ports were located at the nozzle entrance and exit 
planes along the axial centerline of the bottom wall. These 
were used to determine the nozzle exit (test section entrance) 
dynamic pressure, which determines the magnitude of the test 
section pressure gradient (see Appendix A). 

The second nozzle was made from steel and was used witl 
an adverse pressure gradient test section configuration. 

This nozzle, shown in Fig. 1, is two feet long, has an exit 
area of forty nine square inches (seven inches square) and 
a contraction ratio of 11.76. Static pressure ports were 
located on the right side of the nozzle at the inlet and 
exit for the purpose of determining the exit dynamic pressure. 
The wake generator was not mounted in the nozzle because of 
the difficulty encountered in trying to center and align it 
in the metal; rather in this case it was mounted just at 
the entrance of the test section. 
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C. WAKE GENERATORS 



The wake generators were smooth Plexiglas cylinders 
spanning the width of the tunnel. A one inch diameter 
cylinder was used with the steel nozzle and a two inch 
diameter cylinder was used with the Plexiglas nozzle; thus 
"equivalent” turbulent wakes (see Results) were generated 
for each test section configuration. The wake velocity 
defect and turbulence resulted from flow separation of the 
turbulent boundary layers on the cylinders. 

D. TEST SECTION 

As shown in Fig. 3> the test section consisted of two 
parallel , vertical side walls separated by flexible and 
adjustable top and bottom walls, all made of Plexiglas. The 
flexible walls were made from one-eighth inch Plexiglas and 
fastened to the nozzle with piano hinges. Eight adjustmen" 
points, each with two threaded rods, allowed these walls to 
be positioned to any desired area distribution and horizon 
tally levelled so that the test section geometry would re- 
main rectangular throughout. Appendix A contains the devel- 
opment of the relationship between the pressure gradient 
and area distribution, based on a one-dimensional analysis. 
Table I contains the dimensions necessary to set the area 
profiles and form the two pressure distributions used in 
this study. 

The right hand vertical wall of the test section contained 
thirteen slots through which the pitot-static probes were 
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inserted into the flow. These slots were nineteen inches 
high, three-eighths inch wide and spaced five inches apart. 

A teflon sliding seal prevented the influx of air through 
the slot in which the probe was being traversed. A filler 
plug for each unused slot was carefully hand milled and 
polished to make it conform perfectly with the flow side of 
the wall, thus eliminating the possibility of the plugs 
protruding into the flow and inducing separation. A thin 
aluminum strip wedged over each filler plug insured a tight 
fit. 

The left hand vertical wall of the test section contained 
three rectangular access holes equally spaced along the 
axial centerline. These enabled the operator to remove and 
insert the wake generators, measure the probe position 
accurately, and level the lower adjustable wall. Also, along 
the axial centerline of this wall were thirty-three static 
pressure ports spaced one-and-a-half inches apart. A 
manometer bank connected to these ports gave a visual indi 
cation of the linearity and uniformity of the pressure dis- 
tribution. The manometer also served to indicate any regions 
of separation, as did several rows of tufts which were 
attached to the walls of the adverse pressure gradient test 
section. 

E. DIFFUSER AND MOTOR 

The test section was followed by two seven-inch square, 
four feet long steel sections, a diffuser section and a 
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centrifugal blower driven by an electric motor shown in 
Figure 1|. The diffuser was attached to the blower by a 
rubber buffer section which reduced the amount of bivration 
transmitted from the motor to the tunnel. The centrifugal 
blower, operated at three thousand revolutions per minute, 
pulled air through the tunnel and exhausted it vertically 
through a muffler. The total mass flow through the tunnel 
was controlled by a damper mechanism just prior to the 
muffler. 

E. DATA SENSORS 

All data were taken with two calibrated Prandtl-type 
pitot-static probes, one measuring three-and-a-half inches 
and the other eighty-three hundredths inch, stem to tip, 
thus allowing acquisition of data between test section data 
stations. Two U-tube manometers we re connected to the probes 
so that direct readout of local static and dynamic pressure 
was obtained. These manometers, shown in Fig. 5 , could be 
inclined at any angle between two and forty-five degrees, 
providing a wide range of sensitivity. A similar manometer 
was connected to the nozzle static pressure ports so that the 
test section inlet dynamic pressure could be continuously 
monitored. 

The traversing mechanism used to position the pitot- 
static probes in the wake was hand-operated and completely 
detachable from the tunnel. The probe was clamped firmly in 
its movable support fixture on the traverse which was moved 
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vertically by a screw drive mechanism. Probe position could 
be determined to within a hundredth of an inch with the aid 
of a ruled indicator on the traverse. 
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III. EXPERIMENTAL PROCEDURE 



A. TEST SECTION PREPARATION 

The test section was configured to produce both a 
favorable-only pressure gradient (Figure 6), and an adverse- 
favorable pressure gradient (Figure 3) over its entire 
length. For ease in setting the flexible walls to the de- 
sired profile, two reference lines were scribed at the top 
and bottom of the right hand vertical wall of the test sec- 
tion. All dimensions given in Table I are referenced to 
these lines and indicate the distance in inches from these 
reference lines to the flow side of the flexible walls. 

The following procedure was used in setting the flexible 
walls. Proceeding from the test section inlet, the right 
side of each adjustment fixture was set to the desired posi- 
tion and tightened. After all eight fixtures had been set, 
each was re-checked to insure that its setting had not been 
disturbed. Then, with the aid of a bubble level, the left 
side of each adjustment fixture was moved until the wall was 
level. Even though the walls could only be set precisely 
to the desired dimensions at the eight adjustment points, 
the resulting area distribution in the test section never 
varied more than five percent from the desired area. The 
only variance occured when the axial area gradient was large, 
such as at the beginning of the favorable pressure gradient 
profile or in the transition section of the adverse-favorable 
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pressure gradient profile; otherwise the wall profile 
conformed exactly to the desired area distribution. 

B. DATA ACQUISITION 

All data were recorded from the water-filled U-tube 
manometers (Figure 5) connected to the pitot-static tubes 
and nozzles. A scale was placed between the tubes of each 
manometer providing a minimum reading of a hundredth of an 
inch deflection on the slant height. To eliminate any 
possible parallax and to assure consistent readings from 
the manometers, a portable plastic sight unit was placed 
over the manometer scale. This unit restricted one’s line 
of sight to a plane perpendicular to the manometer. 

At each position of the probe in the wake, the manometer 
readings were recorded and plotted on a graph so that it was 
immediately obvious when more data was needed in a particular 
portion of the wake. This technique was particularly advan- 
tageous in determining the exact center of the wake and the 
point where the freestream was encountered. 

The high turbulence intensity in the central core of the 
wake caused large fluctuations in the manometer water columns. 
These fluctuations were often as great as two-tenths inch 
even at manometer inclinations of thirty degrees. By ob- 
serving the fluctuations for several minutes at each probe 
position in the wake, it was possible to ’’average” these 
fluctuations by eye and obtain a reading that was accurate 
to within an estimated five-hundredths of an inch. 
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The inclination (set with an inclinometer) of each of 
the three U-tube manometers was varied according to tunnel 
operating conditions and axial location of the pitot-static 
probe. For a given axial data station and local flow veloc- 
ity, an inclination was chosen that would yield approximately 
a ten to eleven inch difference in water column slant height 
when the probe was positioned in the freestream. This 
assured adequate sensitivity and readability of the water 
columns at the lower velocities of the wake centerline. 

Traversing of the wake by the pitot-static probe was 
always done from the upper free stream edge to the lower 
free stream edge. The free stream was located by moving 
the probe away from the upper wall until the maximum slant 
height difference on the dynamic pressure manometer remained 
unchanged for a few tenths of an inch below the point at 
which it first occurred. This point was taken as the free 
stream. The lower free stream was found in the same way. 
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IV. DISCUSSION OP RESULTS 



A. TEST SECTION 

A qualitative indication of the linearity of the test 
section pressure gradient in both the favorable and the 
adverse-favorable configurations can be seen in Figures 7, 

8. The limitations of the test section in generating an 
exactly linear pressure distribution result from the lack 
of completely flexible upper and lower test section wall 
and from the limited number of wall adjustment points from 
which to set the area profile. Considering these physical 
limitations and the fact that the desired area distribution 
in the test section was derived from one-dimensional analysis 
the resulting actual pressure distribution in the test 
section is quite good. A comprehensive program to calibrate 
the test section was not undertaken since only a qualitative 
idea of its performance was necessary for this investigation. 

B. EXPERIMENT 

The experiment was designed to compare the decay of tw r o 
similar wakes which have had different flow histories. This 
was done, within the physical constraints of the test section 
by using a two-inch diameter wake generator in a favorable 
pressure gradient channel sixty eight inches long and a one- 
inch diameter wake generator in an adverse-favorable pressure 
gradient channel with the favorable pressure gradient 



occupying the downstream thirty four inches (Figures 3 and 
8). The two wakes are compared within the favorable 
pressure gradient regions downstream of the point where 
both wakes have the same velocity defect. The approximate 
solution of the integral momentum equation in Appendix C 
illustrates the possibility that an adverse pressure gradient 
can sustain a wake deficit, such as the one studied here. 

In addition, the wakes are compared at the same relative non- 
dimensional distance into the favorable pressure gradient. 
This was possible because the favorable pressure gradient 
lengths were in the same ratio as the initial wake widths. 

It is apparent from the experimental data in Figure 9 that 
the velocity defect decays rapidly when the flow enters the 
favorable pressure gradient directly from the nozzle. How- 
ever, the velocity defect persists [Figure 10] and then 
decays at a much slower rate [Figure 11] when the flow has 
previously encountered an adverse pressure gradient. Sine 
the rate at which the velocity defect diminishes is direct _y 
related to the eddy viscosity, it can be at least qualita- 
tively concluded from Figures 12,13 that the eddy viscosity 
of a turbulent wake is a function of the flow history and 
its effect on the turbulent field. Classically,, it is 
assumed that the eddy viscosity is purely a function of the 
local mean flow properties, we must conclude that care must 
be exercised when applying classical eddy viscosity models 
to pressure gradient situations. In order to test this 
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conclusion, a numerical analysis of the experimental data 
was conducted. 

C. COMPUTER PROGRAM 

The data analysis computer program uses a forward 
difference scheme to solve the equations of motion for the 
turbulent shear coefficients (CTAU) and the eddy viscosi- 
ties (EDDY) as described in Appendices B and D. The exper- 
imental data obtained from the wind tunnel were used to 
test the program and assure correct processing with no 
coding errors. The numbers generated for CTAU and EDDY were 
unreasonably large and therefore the results are inconclusive. 
It is felt that the increments in the data in the axial 
direction used in finite differencing the data were too 
large for an accurate approximation of the partial deriva- 
tives which appear in the equations (see Appendix B). A 
better method would be to smooth the data and use a central 
difference numerical scheme to approximate the partial 
derivatives. The data taken in the experimental phase of 
this investigation were not adaptable to the central 
difference scheme. Given experimental data with sufficiently 
small increments in the axial direction, the program will 
provide the correct values for the parameters listed in 
Appendix D. 
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V. CONCLUSIONS 



The experiment described in this report was designed to 
test the validity of classical eddy viscosity models in 
free turbulent mixing situations with axial pressure gra- 
dients. The results indicate that those classical models 
which are explicitly independent of the local turbulence 
structure will fail to predict the experimentally observed 
phenomena which are a strong function of the turbulent field. 
In these pressure gradient situations, the local turbulent 
kinetic energy is not a unique function of the local non- 
turbulent mean flow and the influence of the turbulence 
on the form of the Reynolds stresses must be an explicit 
one in any eddy viscosity model. 
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REFERENCED TO SCRIBED LINES 12-INCHES FROM THE TEST SECTION CENTERLINE AT TOP AND BOTTOM OF 
TEST SECTION AND MEASURED TOWARDS TEST SECTION CENTERLINE TO FLOW SIDE OF ADJUSTABLE WALL 




Figure 1 

View of tunnel showing nozzles used with (1) fav- 
orable pressure gradient and (2) adverse-favorable 
pressure gradient test section configurations. 
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Figure 2 

View of Plexiglass nozzle and test section showing 
(1) 2-inch wake generator installed in the nozzle 
and (2) the pitot tube traverse mechanism. 
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Figure 3 

View of Plexiglass test section showing (1) the 
adjustable walls in the adverse-favorable pressure 
gradient configuration, (2) the adjusting rods for 
these walls and (3) the slots through which the pitot 
tubes are inserted. 
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Figure 4 

Overall view of the wind tunnel showing the diffuser, 
centrifugal blower and its muffler at extreme right 
in picture. 
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Figure 5 

View of inclined U-tube manometers used to obtain 
pressure readings in the nozzle and test section. 
The sight unit (1) used to eliminate parallax in 
reading the manometers is shown on the left. 
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Figure 6 

Test section configured for the favorable pressure 
gradient. The milled slot through which the pitot 
tube is inserted can be seen at station *1 . 
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Figure 7 

View of manometer bank connected to the static 
pressure ports along the test section centerline. 
Shown is the linear pressure distribution in the 
test section resulting from the configuration in 
Fig. 6. The horizontal water columns are not 
connected to the tunnel. 
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Figure 8 

View of manometer bank connected to the static 
pressure ports along the test section centerline. 
Shown is the pressure distribution in the test 
section resulting from the configuration in Fig. 3. 
The horizontal water columns are not connected to 
the tunnel. 



36 



STATION NO. 



CO 

Q. 




37 



VELOCITY PROFILES IN A FAVORABLE-ONLY PRESSURE GRADIENT 
(2" DIAMETER WAKE GENERATOR, 68 j " GRADIENT LENGTH) 
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VELOCITY PROFILES IN ADVERSE PRESSURE GRADIENT 
0 " DIAMETER WAKE GENERATOR, 29 \ " GRADIENT 
LENGTH) 
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VELOCITY PROFILES IN A FAVORABLE PRESSURE GRADIENT, FOLLOWING 

AN ADVERSE PRESSURE GRADIENT 

0" WAKE GENERATOR, 29 -^" F. P. G. LENGTH) 
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FIG. 13. COMPARISON OF VELOCITY RECOVERY RATES IN THE TWO LINEAR FAVORABLE 

PRESSURE GRADIENTS 



O FAVORABLE GRADIENT-HIGH SPEED 
O FAVORABLE GRADIENT-LOW SPEED 
A ADVERSE GRADIENT -HIGH SPEED 
□ ADVERSE GRADIENT -LOW SPEED 
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FIG. 12 VELOCITY RECOVERY RATES IN BOTH TEST SECTION CONFIGURATIONS 



APPENDIX A 



DESIGN OF THE TEST SECTION AREA DISTRIBUTION 

The local cross-sectional area of the test section and 
the static pressure of the fluid flowing through the section 
are related through the equation of continuity and Bernoulli' 
equation. For the steady, isentropic flow of an incompres- 
sible fluid, these equations can be written: 

VA = constant (1) 

p + |-pV 2 = constant (2) 

In the log-derivative plane (1) and (2) become 

dV = dA n N 

V “ A 

dp = -pv 2 («•> 

Substitution of (3) into (4) yields 

- , r 2 dA * / c\ 

dp = pV — (5) 

From (1) V = C^/A, where C 1 is the volumetric flow rate at 
the inlet to the test section. Thus, (5) can be written 



( 6 ) 



d£ = dA 

dx p a 3 dx 

After specifying a particular pressure gradient, (6) can 
be solved to determine the required area distribution. 

For a constant pressure gradient, ^ = _c 2 , and inte- 
gration of (6) for this simple case is straightforward. 
The result can be written 



A(x) = 1 

A 1 1 + 6x 



(7) 



where B = C^/q^, A^ is the test section inlet area and q^ 
is the dynamic pressure at the test section inlet. For a 
specified inlet and exit area 




where L is the length of the pressure gradient, and A^ is 
the channel area at x = L. 

The area profiles for the favorable and adverse-favorable 
pressure gradients used during this investigation are shown 
in Figs. 3 and 6, respectively. Table I contains the specific 
dimensions used to set the flexible walls of the test section 
to the desired area profile. Shown for comparison in this 
table are measurements of the upper and lower flexible walls 
as they were set prior to taking data. The differences in 



^3 



dimensions, where they exist, between the theoretical and 
actual wall settings show the degree of conformity of the 
flexible walls to the desired profile. 
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APPENDIX B 



DEVELOPMENT OF AN EQUATION FOR CALCULATION 
OF THE LOCAL TURBULENT SHEAR COEFFICIENT 

The momentum equation governing the mean flow properties 
of a steady, two-dimensional, Incompressible turbulent 
shear flow can be written, using Reynolds 1 averaging 
technique [Ref. 3]: 



u 



3u 

9x 



+ 



3u 1 dp + 
3y “ ” p dx 



i a_r_ 
p 3y 



( 1 ) 



2 

where t = -pu'v* and pu' has been neglected with respect 
to p. The same arguments which Prandtl used to derive 
equation (1) from the Navier-Stokes equations, in the case 
of a boundary layer, are directly applicable to a two- 
dimensional wake in approximately parallel flow through a 
rectangular duct.- Introducing the integrated form of the 
continuity equation 



v 



/<!£> ** 



( 2 ) 



into equation (1) transforms it to 



u 



3u 

9x 



3u / 3u% 

3y o 3x 



dy 



1 dp + 1 3 t 
p dx p 3y 



(3) 



By non-dimensionalizing equation (3) with respect to the 
local freestream velocity, u 0 , and utilizing the product rule 
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of differential calculus, equation (3) becomes 



u 



y rr . 



U 



^ ( n \ < / u \ 2 x /U\ /• / u \ j , u x / x \ ^ / 

u „ u „ x % uT " ( % ) y o % x y " % y ^ 0 ir )d y 

e 



e e 



1 TT 



PU 



_ P + -2L 
2 x 2 



(4) 



2 

where = 2x/pu e and subscripts x and y indicate partial 
differentiation with respect to x and y. Equation (4) is 
now integrated from y = 0 to y = y , an arbitrary point in 
the wake, and relax the condition that equation (*l) be 
satisfied at every point in the flow, allowing it to be 
satisfied on the average between y = 0 and y = y . After 
the necessary integration by parts and considerable manipu- 
lation, the final equation for the turbulent shear coeffi- 
cient, C , can be written 

T 



C T (y T ) 



y T r ,u, 

i L \ ) u X 

0 e e e 






dy 



P U o 



(5) 



Equation (5) was solved numerically with the computer 
program listed in Appendix D. 

The validity of equation (5) can be checked by letting 



u( y T ) - u e ’ "■ 0 

dP _ nu 
dx “ pu e dx 
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Thus C T (y T ) = 0 and equation (5) reduces to 



d_0 

dx 



^e 

dx 



■(20 + 6 *) = 0 



( 6 ) 



where 



• E t 11 ' t’ dy 



6 * 



/ (1 - T~) dy 
0 e 



Equation (6) is the well known momentum integral equation 
in boundary layer theory. 
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APPENDIX C 



SOLUTION TO THE MOMENTUM 
INTEGRAL EQUATION 

The momentum integral equation obtained in Appendix B 
can be integrated in the special case of a constant pressure 
gradient. Prom Eq. (6) in Appendix B 

J n 1 dU 

i£ + ;r -dF (2e + fi *) = 0 (1) 

e 

where 0, 6* are the momentum and displacement thickness 
defined in Appendix B. 

du 

An expression for can be obtained from the inviscid 

Bernoulli equation 




1 dp 
pu e dx 



( 2 ) 



Substitution of Eq. (2) into Eq. (1) yields a differential 
equation in terms of 0, and H = 6*/0 (the shape factor). 



|r(ln 6) - (jjf) (H + 2) = 0 (3) 

PU 



The solution of Eq. (3) is 



= exp[ (fe) (H + 2)dx] (4) 

6 1 0 pu e 2 dx 
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From Appendix A, -£■ = -Bq^ and the integrand of Eq. (4) 
can be written 

3 q i 

- ( 2^ } (H + 2) (5) 

From the results of Appendix A, the ratio is a known function 
of x 



^e 



1 

1 + Bx 



( 6 ) 



After substituting Eqs. (5), (6) into Eq. (4) and integrating, 
the final solution becomes 



1 + 



exp[ / 



BH 



bTT 



dx] (7) 



Thus, experimentally derived values for the shape factor 
can be substituted into Eq. (7), and the resulting value of 
6 can be compared with the experimentally obtained values. 

For the restrictive case, H = constant, we can extend 
this solution to 



6 1 = (1 + 6 x ) 1+H/2 

Further, for H =* 1.0 (a close approximation under the 
conditions of this investigation), Eq. (7) reduces to 
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1 



( 9 ) 



6 1 (1 + fix) 3/2 

This result compares favorable with experimental data and 
serves as a guide for calculating the growth or decay of 
wakes in pressure gradients. 
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APPENDIX D 



DATA REDUCTION PROGRAM 



A . GENERAL 

The following FORTRAN IV computer program uses the 
experimental values of dynamic and static pressure through 
the wake in inches of water and calculates the wake momentum 
thickness, displacement thickness, turbulent shear coeffi- 
cient, eddy viscosity, and Prandtl mixing length. The 
results are printed out as well as plotted on the calcomp 
plotter. 

B. NOMENCLATURE 

The following parameters and arrays are used in the 
program, with their respective dimensions shown in 
parentheses : 

ALV - array used to store interpolated values of 
dynamic pressures (feet) 

AMIXL - local Prandtl mixing length in the wake (feet) 

CTAU - local turbulent shear coefficient for each 
point in the wake 

DELRAT - ratio of succeeding DELSTR’s to that at the 
first data station 

DELSTR - wake displacement thickness (feet) 

DELX - axial distance between succeeding probe tip 
positions (feet) 

DIAMTR - diameter of wake generator (inches) 

EDDY - eddy viscosity 
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EPSIL 



- residual in momentum integral equation when 
experimental values for THETA and SHAPE are 
substituted into it 

LABEL - label for subroutine DRAW 

MN - number of data points laterally through the 
wake 

MX - number of axial data positions 

MY - number of points (always odd) into which the 

wake is divided 

P - array of local static pressures for each point 

in the wake (psf) 

PATM - atmospheric pressure (inches of mercury) 

Q - freestream dynamic pressure (psf) 

SHAPE - wake shape factor (DELSTR/THETA) 

THETA - wake momentum thickness (feet) 

THETAP - inclination of static pressure manometer (degrees) 

THETAV - inclination of dynamic pressure manometer 
(degrees ) 

THRAT - ratio of succeeding THETA 1 s to that at the 
first data station 

TITLE1 - title for subroutine DRAW 

TITLED - title for subroutine DRAW 

U - array of velocities for each point in the 

wake (feet per second) 

URATIO - ratio of local wake velocity to local free- 
stream velocity 

XIN - array of probe traverse positions through the 
wake (inches) 

XOUT - interpolated probe traverse positions for each 
of the MY points in the wake (inches) 

XPOS - probe tip position measured from downstream 
side of the wake generator (inches) 
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YOUT - interpolated manometer readings (for YQ and 
YPS) for each of the MY points through the 
wake (inches) 

YPS - static pressure manometer slant distance 
height difference (inches) 

YQ - dynamic pressure manometer slant distance 
height difference (inches) 



C. SUBROUTINES 

INTEG - uses Simpson’s rule to numerically integrate 
the function F. It must be supplied with K 
points (always odd) equally spaced, DELX 
units apart. 

PIF 2 - second-order interpolation function. X is the 
point for which interpolation is desired. 

XLIST and FLIST contain the N number of abscissas 
and ordinates, respectively, of the table to 
be interpolated. 

DECIDE - insures that the wake is divided into an odd 
number of points spaced 0.1 inches apart. If 
the wake contains an even number of experimental 
data points (found by subtracting the lower 
traverse position from the upper position and 
multiplying by ten), DECIDE adds a point onto 
the upper edge of the wake and sets its value 
equal to the freestream value. The total 
number of points are printed out. 

FILL - insures that the wake at succeeding downstream 
data positions contains as many data points as 
the one before it. Thus division by zero is 
precluded in the computation scheme used for 
EPSIL and CTAU. If the succeeding data station 
has fewer points in its wake, FILL copies the 
uppermost wake point (the freestream point) 
into additional points as necessary to equal the 
number of points at the current position. 

DRAW - SSP 3 subroutine for the IBM 360/67 at the Naval 
Postgraduate School which uses the CALCOMP 
plotter to graph the program output . 
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D. USAGE 



The program will accept up to twelve data positions in 
the streamv/ise direction. If more than that number are 
required, the dimensions of the appropriate arrays (ALP, 
ALV, U , P , A , B , C , D , EDDY , AMI XL , XOUT , URATIO , CTAU , DELSTR , THETA , 
SHAPE , EPSIL , Q , DELRAT , THRAT , THETAV , THETAP , P ATM , XPOS , DELY , 
NUM) will have to be modified accordingly. The program 
is set up for exactly three complete runs (a run is defined 
as the processing of data taken for a particular test 
section configuration and inlet dynamic pressure). If more 
or fewer runs are necessary, the limits for the KK DO-LOOP 
have to be modified accordingly. There is no restriction 
on the number of data points that can be obtained at each 
axial data position. Subscript I always refers to each 
axial position for which data was taken, and subscript J 
always refers to the points through the wake, at a given 
axial data position. 

The following input data is supplied by the user in the 
order presented. Input format is given in parentheses. 

Two cards: TITLE1 (6 a 8) - first graph title 

Two cards: TITLED (6A8) - second graph title 

One card: LABEL (12A4) - curve label. One number in each 

four columns to designate the data position, 
i.e. 1 2 3 etc. 

For each run the following data is needed: 

One card: run identification (up to 47 characters) 
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One 


card : 


number of data positions (12); wake 
diameter, DIAMTR 'inches , (F10.2) 


generator 


One 


or 


more 


cards : 


THETAV, degrees (8F10.1) 




One 


or 


more 


cards : 


THETAP, degrees (8F10.1) 




One 


or 


more 


cards : 


PATM, inches (8F10.2) 




One 


or 


more 


cards : 


XPOS, inches (8F10.2) 




One 


card : 


data position identification (up to 


^7 characters) 


One 


or 


more 


cards : 


XIN s inches (8F10.3) 




One 


or 


more 


cards : 


YQ, inches (8F10.3) 




One 


or 


more 


cares : 


YPS j inches (8F10.3) 





E. PROGRAM LISTING 

A listing of the program which has been described here 
is included on the following six pages. 
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on 



IMPLICIT REAL *4 (A-H,0-Z) 

RE AL* 8 TITLEM12) 

REAL*8 T I TL E 4 ( 12 ) 

REAL*8 LAB ELI / ' DELS'/ 

REAL*8 LABEL2/' THET'/ 

REAL*8 LABEL3/ ' SHAP'/ 

REAL*4 LAB EL ( 1 2 ) 

DIMENSION ALP(12,125),ALV(12,125),U(12,125),P(12,125) 
DI MENS ION A{ 12, 125 J ,B( 12, 1251 ,C( 12, 12 5) ,D( 12, 12 5) 

DIME NS I ON EDDY (12,125) , AM I X L ( 12 , 12 5 ) , XOUT ( 12, 125) 
DIMENSION UR AT ID ( 12,125) ,CTAU( 12, 125) ,YQ (30) , YPS (30) 
DIMENSION UR (125) ,UP( 125) ,UC ( 125) ,UD( 125) ,F( 125) 
DIMENSION YOUT( 125) ,XI N(30) 

DIMENSION AXI S ( 125) 

DIME NS I ON DELSTR ( 12 ), THETA ( 12) , SHAPE ( 12 ) , EPSILC 12) 
DIMENSION Q( 12) , DEL RAT ( 12) , THR AT ( 12 ) ,THETAV(12) 
DIMENSION THETAP (12),PATM( 12) , XPO S( 12 ) , DE L Y( 12) 
DIMENSION NUM (12) 

COMMON/A/ H 

EQUIVALENCE ( F ( 1 ) , Y OUT ( 1 ) , UR ( 1 ) , UP ( 1 ) , UC ( 1 ) , UD( 1 ) ) 
EQUIVALENCE ( A (1,1) , EDDY (1,1) ) 

EQUIVALENCE ( B ( 1 , 1 ) , AMI XL ( 1 , 1 ) ) 

EQUIVALENCE ( C ( 1 , 1 ) , AX IS ( 1 ) ) 

DATA GAMA/62. 32/, RHO/O. 002378/ 

H=0 . 1E0 
Y =H/12 • 0 



50 FORMAT (' 

1 • ) 

51 FORMAT ( I 2, F 10. 2 ) 

52 FORMAT (8F10.1 ) 

53 FORMAT ( 8F1 0.2 ) 

54 FORMAT ( 1 
1 ' ) 

55 FORMAT (12) 

56 FORMAT ( 8 F10 . 3) 

57 FORMAT (/////// 1 *,T20,6A8) 

58 FORMAT! /////// ' *,T20,12A4) 

70 FORMAT (//////////' ' , T 33 , ' ECHO-CHEC K DATA'///) 

71 FORMAT ( / / ' ' ,T20,' THETAV' , (' ' , T33 , 8 F 1 J . 2/ / ) ) 

72 FORMAT!//' ', T20 ,' THETAP ',( ' ' , T33 , 8F 1 0. 2/ / ) ) 

73 FORMAT ( / / ' ' , T 20 , ' P ATM ' , ( ' ' , T33 , 8F 10. 2/ / ) ) 

74 FORMAT { // ' ' , T 2 0 , ' X POS * , ( ' ' , T33 , 8F1 0 . 2/ / ) ) 

75 FORMAT ( /// ' ' , T 53 , ' MN = ' , I 3/ / / ) 

76 FORMAT ( // ' ' ,T 1 0 , ' X IN' , ( ' ' , T 2 1 , 8 F 1 0. 3/ / ) ) 

77 FORMAT ( / / ' ' , T 1 0 , ' Y Q ' , ( ' ' ,T2 1 , 8F1 0 .3/ / ) ) 

78 FORM AT ( / / ' ' , T 1 0 , ' YPS ' , ( ' ' , T2 1 , 8F 1 0. 3/ / ) ) 

79 FORMAT ( ' 1' ,T50, '***RESULTS OF DATA REDUCTION***') 

80 FORMAT ( //////' ',T50, 'AXIAL STA T I ON ' , I 2 , ' OUT OF ' , I 2 , 

l'ST AT IONS'//// / ) 

81 FORMAT ( ' ' ,T18 , ' DEL ST AR=' , G12 .5, 5X , 'THET A= ' ,G12. 5, 5X, 

1 *SHAPE=' ,G12.5, 5X, 'EPSILON^' ,G12.5) 

82 FORMAT (//////' ' ,T8 ,' POINT ', T30 , *U ', T48, ' U/UE ', T70, « P ' 

1 , T 8 8 , ' C T A U ' ) 

83 FORMAT ( '0',T9,I3,T24,G12.5,T44,G12.5,T64,G12.5,T84, 
1G12.5) 

84 FORMAT ( ' 1 ' , T 44 ,'*** SUMMARY OF INTEGRATED WAKE PROPER', 
l'TIES***' ) 

8 5 FORMAT (///// ' ' , T3 , ' AX I AL ' , T 7 5 , 'THETA'/' ' ,T2, 'STATION 

1« ,T16, ' X/D' ,T29, 'DEL STAR' , T44, • DELRATIO' ,T60 ,' THETA' , 
2T75, 'RATIO ' ,T90, ' SHAPE ' ,T105, 'EPSILON') 

86 FORMAT ('O' ,T4,I3,T11,G12.5,T26,G12.5,T41,G12.5,T56, 
1G12.5* T71.G12. 5, T86 ,G12. 5 , T102 ,G12. 5) 

87 FORMAT {' 1' ,T49, '***REYNOLD"S MIXING LENGTH DATA***') 

88 FORMAT (//////' ' ,T6 1 , ' ST AT ION' , 13/// ' ' , T 44 , ' PO I NT ' , 

IT 59 , 'EDDY V I SC 0 S I T Y ' , T 79 , ' MI X I NG LENGTH') 

89 FORMAT (' ' ,T45, I3,T60,G12.5,T80,G12.5) 

90 FORMAT ( 6A8 ) 

91 FORMAT { / / ' ' , ' L AST= ' , I 3 ) 

92 FORMAT (12A4) 



56 



non non o no 



READ ( 5 , 90 J TITLE1 
RE AD ( 5 , 90 ) TITLE4 
READ (5*92) LABEL 

ZERO ALL ARRAYS BEFORE COMPUTATION BEGINS 
F ( J ) =Z 



DO 1000 KK = 1 

1 = 0.0 

DO 1 J = 1 , 125 
DO 1 1=1,12 
XOUT ( I , J)=Z 
ALP ( I, J)=Z 
ALV ( I, J)=Z 
U ( I , J)=Z 
URAT I 0 ( I ,J) = 
P( I,J) = Z 
A(I , J)=Z 
B( I , JJ =Z 
C( I, J ) = Z 
D( I, J)=Z 
CTAU( I , J)=Z 
DEL STR ( I )= Z 
THETA! I ) = Z 
SHAPE ( I ) =Z 
EPS I L ( I J = Z 
0 ( I ) =Z 
DE LR AT ( I )=Z 
THRAT! I )=Z 
THETAV! I ) = Z 
THETAP! I ) =Z 
PATM! I )=Z 
XPOS! I )=Z 
DEL Y( I ) =Z 
CONT INUE 
DO 2 1=1 ,30 
X I N ! I)=Z 

y o ( ij = z 

YPS ( I ) = Z 
CONTINUE 
READ ( 5 , 50) 

RE AD ( 5 , 5 1 ) 
READ ( 5 , 52) 

R EAD ( 5,52) 
READ! 5,53) 
READ( 5, 53) 

W R I T E ( 6 , 5 0 ) 
WRITE! 6,70) 
WRITE(6,71) 
WRITE (6 ,72) 
WRI TE ( 6 , 73 ) 
WRITE! 6,74) 
WRITE !6, 57) 
WRITE! 6, 57) 
WRITE (6,58) 



,3 



MX, D I AMTR 
(THETAV! I) ,1=1 , MX ) 
(THETAP ( I ) , 1 = 1, MX) 
( PATM! I ) ,1=1 , MX ) 

( XPOS! I ) ,1 =1,MX) 



( THETAV! I ) ,1 = 1, MX) 
(THETAP! I ) , 1 = 1, MX ) 
( PATM! I) ,1=1, MX) 

( XPOS! I ) , 1 = 1, MX) 

( TITLE1! I ) , 1=1 ,6) 

( TITLE4! I ) ,1=1,6) 

( LABEL! I ) , 1=1, 12) 



BEGIN COMPUTATION FOR EACH STATION 

MX1 = MX— 1 
DO 6 1 = 1, MX 

THETAV ( I )=THET AV ( I )/57 .296 
THETAP! I ) =THET AP ( I )/5 7.296 
PATM! I ) = P A T M ( I )*70.73 
XPOS! I)=XPOS(I ) / 1 2 . 0 

READ IN MANOMETER HEIGHTS 



READ! 5,54) 
READ! 5, 55) 
R EAD ( 5 , 56 ) 
READ! 5 , 56) 
READ(5, 56) 



m 

(X IN! J ) , J = 1 , MN ) 
( Y 0 ( J ) , J = 1 , MN) 
( YPS! J ) , J = 1,MN) 
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non non non non non ono no 



ECHO-CHECK INPUT 

WRITE (6,54) 

WRITE(6,75) MN 

WRITE (6, 76) (XIN(J),J=1,MN) 

WRITE (6, 77) (YOU ) , J=1,MN) 

WRITE (6,78) (YPS(J),J=1,MN) 

CALL DECIDE(XIN,MN,MY) 

NUM ( I ) = MY 

INTERPOLATE FOR MISSING DATA POINTS INSIDE WAKE 

XOUT ( I » 1 ) =XI N( 1) 

DO 21 J = 2» MY 

XOUT (I , J)=XOUT( 1,1 )+FLOAT (J— 1 )#H 
21 CONTINUE 

DO 3 J = 1 » MY 
X=XOUT ( I ,J) 

YOUT (J)-PIF2(X,XIN,MN,YQ) 

ALV ( I » J ) =Y OUT ( J ) / 1 2 .0 

U( I , J ) = SORT ( 2. O^GAMA^SI N(THETAV(I ) )*ALV ( I , JJ/RHO) 

3 CONTINUE 

Q ( I ) =0 • S^RHO^U ( I ,MY )**2 

CALCULATE DISPLACEMENT THICKNESS 
DO 31 J= 1 , MY 

URAT I 0 ( I , J)=U( I , J)/U( I,MY) 

F ( J ) =1 • 0—URAT I 0 ( I , J ) 

31 CONTINUE 

CALL INTEG(F, MY, Y, RESULT) 

DE LSTR { I ) = RESULT 

DELRAT ( I ) = DELS TR ( I ) /DE L STR ( 1 ) 

CALCULATE MOMENTUM THICKNESS 

DO 4 J = 1 , MY 

F(J)=URATIO(I , J)*( 1.0-URATI 0( I , J) ) 

4 CONTINUE 

CALL INTEG(F,MY,Y, RESULT) 

THETA ( I ) =RE SUL T 

THRAT ( 1 ) =THETA ( I ) /THETA ( 1) 

CALCULATE SHAPE FACTOR IN MOMENTUM INTEGRAL 

SHAPEt I ) =DELST R ( I J/THET A( I ) 

INTERPOLATE P-STATIC MANOMETER READINGS 

DO 5 J=1 t MY 
X=XOUT ( I , J ) 

Y OUT ( J )=PIF2 (X,XIN,MN,-YPS J 
ALP (I , J)=YOUT( JJ/12.0 

P( I » J ) = P ATM ( I ) - GAMA*SIN(THETAP( I ) )*ALP( I , J) 

5 CONTINUE 

6 CONTINUE 

CALCULATE RESIDUAL IN MOMENTUM INTEGRAL EQUATION 

DO 7 I = 1 , M X 1 

DE LX = X POS ( 1 + 1 ) -X POS { I ) 

M Y=NUM ( I ) 

IF(NUM( I ).GT.NUM(I+1) ) CALL FI LL ( NUM, I , U ) 

El=( THETA ( 1+1) -THETAt I ) )/DELX 
E2 = THETA( I + 1)*( SHAPEt I + D+2.0) 

E3=(U( 1+1,MY)-U( I, MY) ) /(DELX*U( I+1,MY) ) 

EPSIL(I + 1) = El + E2* E3 

7 CONTINUE 

EPS IL ( 1 ) =0 • 0 
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non non nn 



13 

14 



15 



16 

17 



CALCULATE LOCAL TURBULENT SHEAR COEFFICIENT IN WAKE 

DO 10 1=1, MX1 

DE LX = X POS ( 1 + 1) -XPGS (I) 

M Y=NUM ( 1 + 1 ) 

IF(NUM( I ) .GT.NUM ( 1+1) ) CALL F I L L ( NUM , I , A L V ) 

DO 10 K=1 , MY 
DO 9 J=1,K 

A ( 1+1, J )=2.0*JRATI0( 1+1, J )-URATI0( I+1,K) 

B ( I+1,J)=(URATI0( I + 1 * J )-URATI0( I , J ) )/ DELX 
C 1 = UR AT 1 0 ( 1+1, J) /U( 1+1 ,MY) 

C2= ( U ( 1+1, MY ) -U ( I, MY ))/ DELX 
C(I+1, J)=C1*C2 

Dl= S IN ( THE TAP ( 1 + 1) ) / S I N ( THE TAV ( I + 1 ) ) 

D2=( ALP ( I , J )— ALP { 1 + 1, J ) ) / (DELX*ALV( I + 1,MY) ) 

D(I+1, J)=D1*D2 

F( J)=A( 1+1, J)*(B(I+1,J)+C( I+1,J))+DC 1+1, JJ/2.0 
9 CONTINUE 

CALL INTEG(F,K,Y, RESULT) 

CT A J ( I+1,K )=2*0*RESULT 
10 CONTINUE 
WRITE! 6, 79) 

DO 14 1=1, MX 
WRITE! 6, 80) I, MX 

WRITE! 6, 81 ) DEL STR ( I ), THETA! I ) , SHAPE! I) , EPSIL (I ) 

WRITE (6,82 ) 

MY-NUM! I ) 

DO 13 J = 1, MY 

WRITE! 6, 83) J , U ( I , J ) , UR AT 1 0 ( I , J ) , P ( I,J) ,CTAU( I, J) 

CONTINUE 

CONT INUE 

WRI TE ( 6 , 84 ) 

WRITE(6, 85) 

DO 15 1=1, MX 

XD = 1 2 . 0*XPOS C I ) / D I AMTR 

WRITE! 6, 86) I , XO , D E L STR ( I ) ,DELRAT ! I ) , TH ET A ( I ) , THR AT ( I ) 
1, SHAPE! I), EPSIL! I) 

CONTINUE 

CALCULATE EDDY VISCOSITY AND REYNOLD'S MIXING LENGTH 

WRITE(6,87) 

DO 17 I =2, MX 

WR I TE ( 6 , 88 ) I 

M I N= ( (NUM! I )-l )/2)+10 

MAX = NU M ( I ) — 1 0 

DO 16 J = MI N , MA X 

E 1 = ABS ( CT AU ( I, J)*Q( I) ) 

E 2 = ABS ( U ( I , J + l )-U( I ,J-1 ) ) 

E3=2. 0"H 

EDDY ( I,J)=(E1/RH0)/(E2/E3) 

AMI XL! I ,J)=SQRT(E1/(RH0*(E2/E3)**2) ) 

WRITE! 6 ,89 ) J , EDD Y ( I , J ) , AMI XL { I , J ) 

CCNTINUE 
CONTINUE 

PLOT PACKAGE BEGINS — MAIN PROGRAM ENDS. 



A A = 0 . 

BB=0 • 

IF (KK.EQ.3) A A = 2 . OE 00 
IF (KK.EQ.2) B B= 1 * OEO 1 
DO 100 1=1, MX 
100 XPOS! I ) =12.^ XPOS! I J/DIAMTR 

CALL DRAW! MX, XPOS, SHAPE, 1,5, LAB EL3, TITLE 1,0., AA, 0,0,0, 

10,9,6,0, LAST) 

CALL DRAW! MX , X PO S, THR AT , 2 , 4 , LAB E L2 ,TITLE2, 0,0, 0,0,0, 

10, 9, 6,0, LAST) 

CALL DRAW! MX , X POS , D E LR AT , 3 ,3, LABE LI ,TITLE1, 0,0, 0,0,0, 

10,9,6,0, LAST) 
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101 

102 

1000 



c 



1 



2 



9 

30 

10 
11 

12 

13 

14 

15 

16 

17 

18 



WRITE (6 ,91 ) LAST 
DO 102 1=1, MX 
MY=NUM ( I ) 

IF(I.EO.l) K=1 
I F ( I .GT.l.AND. I . LT.MX) K=2 
IF(I.EO.MX) K= 3 
DO 101 J=1 , MY 
UR ( J ) =U ( I, J J 

AX I S ( J ) =XOUT ( I , J )/ D IAMTR 

CONTINUE 

LABEL4=L AB EL ( I ) 

CALL DRAW ( MY , AX IS , UR, K, 0, LABEL4, T ITLE4, 0. , BB, 0, 0, 0, 0, 
19, 15,0, LAST) 

CONT INUE 

WRITE(6,91 ) LAST 

CONTINUE 

STOP 

END 



SUBROUTINE I NT EG ( F , K,DELX, RESULT) 
IMPLICIT REAL*4 (A-H,0-Z) 
DIMENSION F ( 1 ) 

INTEGRATION BY SIMPSONS RULE 
K 1 = K~1 
A = F ( 1 ) + F ( K ) 

B= 0 . 0 

DO 1 1=2, Kl, 2 
SUM=4. 0*F( I ) 

B=B+SUM 
CONT INUE 
C= 0. 0 

DO 2 1=3, K, 2 
SUM=2.0*F( I ) 

C=C+SUM 

CONTINUE 

RESULT=DELX*(A+B+C)/3.0 

RETURN 

END 



FUNCTION PIF2(X,XLIST,N,FLIST) 

DIMENSION XLIST( 1) ,FLI ST( 1) 

BLIF(P,Q,R,S,T )=(( Q— P ) * ( S-T )/(R-Q)+S) 

IF(X— XLIST(N) ) 2,2,1 
I = N— 1 
GOTO 5 

I F ( X-XL I ST ( 1 ) ) 4,4,6 

1 = 1 
K = 1 

GOTO 30 
K=2 

DO 8 1=1 ,N 

IF(X— XLIST(I)) 9,9,8 
CONT INUE 
I =N 
1 = 1-1 

BLIF1 = BLIF(X,XLIST(I) ,XLIST( 1 + 1) ,FLIST( I ) ,FLIST( 1+1) ) 
IF(K-l) 11,11,12 
P I F2=BL I F 1 
RETURN 

IF ( ( I + 2-N) ) 13,13,16 

I F ( ( I-D-1) 15,14,14 

IF(A3S(XL1ST(I-1 )— X ) - ABS ( XL I ST (I +2 ) -X ) ) 16,15,15 
L = I + 2 
GOTO 17 
L=I-1 

BL I F2 = BL IF ( X , XL I ST (I) , X LI ST ( L ) , F L I ST ( I ) , F L I ST ( L ) ) 
PIF2=BLIF(X,XLIST( I+1),XLIST(L) ,BLIF1,BLIF2) 

RETURN 

END 
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on 



SUBROUTINE DEC I D E ( X , M , N ) 

COMMON/A/ H 
DIMENSION X ( 1 ) 

10 FORMAT ( ///' • ,T10, ' NUMBER OF POINTS IN WAKE=',I4) 

A= 1 • 0/H 

N= I F IX ( { X { M )— X { 1 ) ) * A ) 

L = 1 

DO 5 K=1,N 
GOTO (1,2), L 

1 L = 2 

5 CONTINUE 

C N IS AN ODD NUMBER OK 

WR I TE ( 6 • 1 0 ) N 
RETURN 

2 L=1 

IF(K-N) 5,3,3 

3 N=N+ 1 

C N WAS EVEN AND WE FORCE IT TO BE ODD 

X( M ) =X ( M)+H 

THE EXPERIMENTAL DATA SHOULD BE NEARLY THE SAME FOR 
THIS POINT AS FOR THE ORIGINAL. 

WR I TE ( 6 , 10 ) N 
RETURN 
END 



SUBROUTINE F I L L ( NUM , I , U ) 

IMPLICIT REAL*4 (A-H,0-Z) 

DIMENSION NUM{ 1) , U ( 12 , 1 ) 

M=NUM ( I ) 

N = NU M ( 1+1) 

DO 1 J=N,M 

U( 1 + 1, J )=U( 1+1 ,N) 

1 CONTINUE 

WR I T E ( 6 , 2 ) M, N 

WRITE (6, 3) ( J ( 1+1, J) , J=N,M) 

2 FORMAT { « O' ,T2J , • M= ' , I 3 , 1 OX , • N= • , 1 3/ / ) 

3 FORMAT ( »0« , 12G10.4/ ) 

RETURN 

END 
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